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Abstract 

We investigated the dependence of the optical depth r of the microlensing events 
on model parameters of the Galactic halo. We only consider Galactic mass models in 
which the rotation curve inside the Sun is compatible with the observation and LMC 
is bound to the Galaxy. It is found that r varies up to a factor 2.5 from the standard 
spherical and flat rotation halo model. This implies that only the most centrally 
concentrated halo model can be consistent with the observation if the halo consists 
of only MACHOs. We also calculate the power x of IMF of MACHO consistent with 
Tyson's CCD survey as well as Bahcall et al. 's observation by HST. It is found that 
x is greater than 5. 

1 Introduction 

Gravitational microlensing events are detected recently by three collaborations: MA- 
CHO (Alcock et al. 1993 [3]) and EROS (Aubourg et al. 1993 [8], Aubourg et al. 1995 
[9]) for LMC events, as well as OGLE (Udalski et al. (1993) [42]) and MACHO (Bennett 
et al. 1994 [13]) for bulge events. From these results we can discuss about the nature of 
the missing mass in our Galaxy (Paczyhski 1986 [30]). 

The key quantity in this problem is the optical depth r of the microlensing. This 
quantity is the instantaneous probability that the event is occurring when we observe 
background stars randomly. Observationally r is derived from the number of the observed 
stars, the mean event duration and the event rate (Paczyhski 1986[30]). For LMC events, 
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MACHO collaboration concluded that r = 8.0tg 4 X 10~ 8 in 68% confidence level(Alcock 
et al. (1995) [5]). On the other hand, several authors (e.g. Paczyhski 1986 [30] Griest 
1991 [24]) estimated r for the halo model with flat rotation curve in the outer region of 
the Galaxy and obtained r ~ 5 X 10~ 7 , which suggests a MACHO fraction / = 0. 21^14 
and all of the missing mass of our Galaxy may not be MACHOs. For bulge events, the 
observational value of r is ~ 3 X 10~ 6 e _1 (Udalski 1994 [43]) by OGLE collaboration, 
where e is the efficiency of the observation, and r = 3.0to [9 X 10~ 6 by MACHO collaboration 
(Bennett 1994 [13]). Paczyhski (1991) [31], Griest et al. (1991) [25] and Kiraga & Paczyhski 
(1994) [26] estimated r for the bulge as 0.1 ~ 1 X 10~ 6 , which is at least factor 3 smaller 
than the observational data, that is, / ~ 3. The fraction, /, however, depends on the 
theoretical estimate of r i.e. models of our Galaxy. So it is important to estimate the 
model dependence of r. The main purpose of this paper is to discuss theoretical r in more 
detail. 

Although the model with the flat rotation curve is frequently taken as the mass distri- 
bution in the halo, at present we can only say that the Galactic rotation curve is essentially 
flat only inside the solar neighborhood (e.g. figure 2 in Fich & Tremaine 1991 [18]). In the 
outer region we have no definite rotation curve at present although we can impose some 
constraints as will be discussed in section 2. Many other spiral galaxies have flat rotation 
curve up to the outer most region. However, it is reported that spiral galaxies with its 
exponential disk scale length less than 3.5kpc have declining rotation curves (Casertano 
& van Gorkom 1991 [15]). For our Galaxy its scale length seems to be marginal, i.e., 
~ 3.5kpc. Our Galaxy may have non-flat rotation curve beyond solar neighborhood. 

The shape of the halo is another point to be considered. It is suggested from N-body 
simulations of the galaxy formation that the halo may be nonspherical (e.g. Aarseth & 
Binney 1978, [1], Aguilar & Merritt 1990 [2], Binney 1994 [14] and references therein). 
Sackett & Gould (1993) [36] and Frieman & Scoccimarro (1994) [19] discussed that the 
ratio of the optical depth toward SMC to LMC is a good probe for the shape of the halo. 
In this paper, we also investigate the dependence of r on the shape. 

As a model of the halo we take the model of Evans (1994) [17] that is a power-law model 
with a rising, flat or falling rotation curve. We impose some constraints in this model and 
calculate the dependence of the optical depth r on model parameters. The same model 
has been used by Alcock et al. (1994) [4]. They concluded that r changes up to a factor 
10. In this paper we use more stringent constraints than theirs, which will be shown later. 

MACHOs may be the low-mass stars. Richer & Fahlman (1992) [33] suggested that the 
IMF of the low-mass star less massive than ~ O.5Af in the Galactic spheroid stars is much 
steeper than the Salpeter's IMF. Bahcall et al. (1994) [12] observed recently halo stars in a 
high-latitude region by HST and concluded that if the dark halo consists of low-mass stars 
they must have mass less than hydrogen-burning limit. Using data of the CCD survey of 
Tyson (1988) [40] and HST observation by Bahcall et al. (1994), we will discuss on the 
constraints to the power of IMF assuming the power law IMF. 

The plan of this paper is as follows. In section 2 we will show the Galactic model and 
impose the constraint on it. In section 3 we will derive the dependence of r on our model 
parameter. In section 4 we discuss the power of IMF consistent with the observation. 
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Section 5 is devoted to discussions. 



2 The Galactic model 



As a spherical or spheroidal Galactic model we take an axisymmetric power-law model 
(Evans 1994 [17]). In this model the gravitational potential is given in the cylindrical 
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,z as 
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where v 0} R C} and q are the normalization of potential, the core radius and the axis ratio 
of equipotential, respectively. The rotation velocity in the equatorial plane is 
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For the other components of the Galaxy we take the model by Bahcall, Schmidt, & 
Soneira (1982) [11] with the exponential disk, the r? spheroid and the central mass concen- 
tration. The surface density of known matter in the Galactic disk at the solar neighborhood 
is H id = 48 ± 8M pc" 2 ([7]). Bahcall, Flynn & Gould (1992) [10] have claimed that the 
total surface density E = 88Af pc~ 2 , implying that there is large amount of the disk DM. 
On the other hand, Kuijken & Gilmore (1989) [28] concluded that E = 46 ± 9Af pc~ 2 , 
and there is no evidence for the disk DM. Since the existence of disk DM is not established, 
we regard E , or the local mass density of the disk p^o = Y< /2zt l} where Zh{= 300pc) is the 
scale height of the disk, as a parameter. We consider two cases, the heavy disk (with disk 
DM) and the light disk (without disk DM), respectively. The density profile of the disk, 
the spheroid and the central mass concentration are respectively 
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where r, i?o(= 8.5kpc), and Rd(= 3.5kpc) are the distance from the Galactic center, the 
galactocentric radius of the sun, and the disk scale length, respectively. As the total mass 
density p we take p = p^ + p<i + p sp h + Pcen- The rotation curve is determined by the sum 
of the gravitational potential of the above four components. 

Now we have five parameters in our Galactic model: R c (the core radius of the halo), 
v (the normalization of potential due to halo), e (the asymptotic ellipticity of the halo), 
f3 (the slope of the rotation curve), and p^o (the local mass density of the disk). We have, 
however, some observational constraints to these parameters: 

1. Although as discussed in section 1 we have no definite rotation curve for R > i? , for 
R < Ro the nearly flat rotation curve is observed. We take this fact into account. We 
fix v c (R = R ) = 225km/s and v c (R = 2kpc) = 200km/s (e.g. [18]) on the Galactic 
plane. This constraint determines the value of v and R c for given values of f3 and 
Pdo- 

2. Kroupa, Roser & Bastian (1994) [27] have observed the proper motion of LMC and 
found that the mean galactocentric space motion vector is (—279 ± 165, —218 ± 
79, +85 ± 122)km/s. Adopting this proper motion, we require that LMC is gravita- 
tionally bound to the Galaxy. This restricts allowed range of /3, e. Too large f3 is not 
possible. 

3. Thomas (1989) [39] observed the motion of the Galactic globular clusters in the halo 
and got the velocity dispersion along the line-of sight direction as 140 ± 37km/s and 
the tangential direction as 223±132km/s for the globular clusters with galactocentric 
radius between ~ 7kpc and ~ 20kpc. Our rotation curve must be consistent with 
this result. 

Under the constraint 1 the rotation curve on the equatorial plane is shown in Figure 1. 
Note that this curve is independent of q since the potential at z=0 is independent of q (see 
equation 2). Using constraints 2 and 3, we have the allowed range of f3 as approximately 
— 1 < ft < 0.6 for the heavy disk and —0.8 < f3 < 0.6 for the light disk. As for e, we 
consider E0~E7 halo (0 < e < 0.7) (Rix 1994 [34]). In reality, the constraint 3 does not 
play an important role. 

The mass of the Galaxy in the above range of parameters is 2 X lO n Af to 2 X 1O 12 M 
within r = 55kpc. To r = lOOkpc it is 8 X lO n Af to 8 X 1O 12 M . This value is consistent 
with the observational estimates (e.g. Ashman 1992 [7]). 

3 The optical depth of microlensing 

In this section we calculate the optical depth of microlensing in our model for LMC 
events and bulge events, respectively. The optical depth of the microlensing r is the 
probability for MACHOs to be in the microlensing tube ([24]) and written as ([24], [30]) 




(9) 
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Figure 1: Rotation curves on the equatorial plane of our Galactic model for specified values 
of f3. (a) is for the heavy disk, while (b) is for the light disk. 
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where Dl is the distance to the lens object from the Sun, m is the mass of the object, Ds 

l/2 

is the distance to the source and Re = [(4Gm/c 2 )(l — Dl/ Ds)] is the Einstein radius. 
We will calculate r for uj = 1. 

In Figures 2 and 3 we show the cumulative contribution f c to r defined by 

MD L ) = I*" P -^^^ T Rl{D L )dDL I f DS P -^^^ T Rl{D L )dDL . (10) 
Jo m j Jo m 

In Figure 2 we show f c for E0 halo. The contribution to r mainly comes from the halo 
component. We also see that the mass distribution in the relatively neighbor region i.e. 
Dl < 20kpc from the sun mainly determines the optical depth. In Figure 3, we show f c for 
E6 halo. The increase of e makes the contributing region more nearer and narrower. We 
show in Figure 4 f c for bulge events with the same notation as for LMC events. For bulge 
events, the contribution is almost independent of (3 and e. The main contribution comes 
from the disk component for the heavy disk, while for the light disk the contribution from 
the disk and from the halo is comparable. 

In Figure 5 we show r to the direction toward LMC with Ds = 55kpc. From this figure 
we see that r does not change more than a factor 2.5 from the model with the spherical 
halo and flat rotation curve (e = 0,/3 = 0), as we change /3, e, and p^o- This result comes 
from the fact that the rotation curve for R < R almost fixes the structure of the halo 
within Dl < 20kpc (or r < 20kpc) where the contribution to r is more than 60 %. We 
also see that r does not depend on e strongly. This is consistent with the results of Sackett 
& Gould (1993) [36] and Frieman & Scoccimarro (1994) [19], although they used a simpler 
model of the Galaxy than ours. 

Using the Galactic model with the halo and disk models as ours, without the spheroid 
and the central component, Alcock et al. (1994) [4] concluded that r changes of an order of 
magnitude. They imposed the constraint that the rotation velocity lies between 180km/s 
and 250km/s at R and 2R . But in their model it is not clear that the inner rotation curve 
is reproduced. Under present constraints in this paper, the rotation curve of the Galaxy 
for r < 15kpc is almost independent of parameters. Since this region has the dominant 
contribution to r as we can see in Figure 2 and 3, the parameter dependence of r is much 
smaller than that of Alcock et al. (1994) 

We also show r to the direction towards Baade's window with Ds = Ro in Figure 6. r is 
almost independent of (3 and e due to the constraints to the rotation curve, r ~ 1.2 X 10~ 6 
for the heavy disk and r ~ 8 X 10~ 7 for the light disk. As mentioned in section 1, MACHO 
collaboration observed that the optical depth for bulge events is 3.0lJ;| X 10~ 6 , which is a 
factor 3 ~ 6 greater than the previous theoretical estimates of Paczyhski (1991) [31] and 
Griest et al. (1991) [25] taking the stellar component of the disk. Our result implies that 
even if we take into account the halo and spheroid component, the optical depth is at least 
a factor 2.5 less than that observed. 

Finally, we show the model dependence of the ratio of the optical depth for SMC and 
LMC, tsmc /tlmc 1 in Figure 7. tsmc I tlmc increases as (3 increases from 0. This is 
because of the following reason. The line of sight toward SMC passes through nearer to 
the Galactic center. For (3 > the outer halo is not important as in the case of (3 = 0. On 
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Figure 2: The cumulative contribution to the optical depth as a function of Dl for LMC 
events for asymptotically EO halo. The disk model is (a)heavy disk (b)light disk 
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Figure 3: The cumulative contribution to the optical depth as a function of Dl for LMC 
events for asymptotically E7 halo. The disk model is (a)heavy disk (b)light disk 
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Figure 4: The cumulative contribution to the optical depth as a function of Dl for bulge 
events for asymptotically EO halo. The disk model is (a)heavy disk (b)light disk 
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Figure 5: The optical depth towards the direction ol LMC as a lunction ol /3, (a) heavy 
disk and (b) light disk. 
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Figure 6: The optical depth of the microlensing towards the direction of Baade's window, 
as a function of f3. (a) heavy disk and (b) light disk. 
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the other hand when f3 < and e > the compression of the less centrally concentrated 
halo along the ^-direction makes the density in the dominant region for the microlensing 
increase specially for SMC, so that tsmc I tlmc increases in this case also. 

From Figure 7 we see for e < 0.3, tsmc I tlmc does not depend on f3 so that it may be a 
tracer of halo flattening, but for e > 0.3 the ratio changes more as a function of f3. So from 
tsmc /tlmc i we obtain a certain combination of flatness and f3. Sackett & Gould (1993) 
[36] and Frieman & Scoccimarro (1994) [19] used only the model with the flat rotation 
curve and concluded that tsmc I tlmc is a good probe of the flattening of the halo if the 
halo-disk tilt angle is ~ 0. Our result shows that even if the tilt angle is 0, we may not 
determine the flattening directly without the knowledge of f3. 

4 The stellar number count and the halo IMF 



Now we calculate stellar number counts in our model and argue on the IMF in the halo. 
Richer & Fahlman (1992) [33] observed low mass stars down to O.14M and suggested that 
the slope of the IMF is steeper than Salpeter's one. More recently Bahcall et al. (1994) 
performed star count of high latitude field by HST and found fewer faint red stars. This 
suggests that the low-mass stars less massive than the hydrogen-burning limit, O.O8Af 
may dominates in the halo, if the dark halo or the part of the dark halo consists of the low- 
mass stars. This is consistent with the statistical mass estimate of MACHOs (Sutherland 
et al. 1994 [38]). 

We now assume ~ 20% of halo consists of MACHO (Alcock et al. 1995 [5]) and ask 
what its IMF would be. We assume power law IMF for MACHO. We use the data of 
Tyson's CCD survey(Tyson 1988 [40], Tyson & Seitzer 1988 [41]). He observed 12 high- 
latitude fields where there are no stars or galaxies brighter than Bj = 20 and R = 19 exist, 
and found about 50 stellar-like objects in a field with 2.6 X 4.6arcmin 2 wide. This number 
count provide a constraint to the halo IMF if the dark halo consists of stars and brown 
dwarfs. 

If the fraction / of the halo consists of MACHO in our model of the Galaxy the expected 
star counts N in Tyson's CCD survey is calculated as 

N = dtt [ D 2 dD I dMp<P x / (11) 

JD>D mt „ J 

= at D^wi AlW fflr« x/ , (12) 

JD>D min Jm<m Um \ dlVl ) Om 

where m is the apparent magnitude of the star, D is the distance to the star, dft is the solid 
angle of the field, D m i n = lOOpc is the minimum distance from which we begin to count 
stars, mu m is the limiting magnitude, and L is the V band luminosity of the star. IMF <f> is 
normalized as / cf>dM = 1. L(M) is the mass-luminosity relation. This relation is derived 
by D'Antona (1987) [16] for the population II stars of M ^ 1M© and approximately given 
as L(M) oc M 2 ' 7 . Now we neglect the contribution of massive stars to N and we assume 
M < O.O8M does not shine. 
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Figure 7: tsmc I tlmc as a function of f3 for various e. (a) is for the heavy disk and (b) is 
for the light disk. 



13 



Figure 8: The total number count N as lunction ol f3 lor various x lor the direction ol 
Tyson's SGP (/ = 312.5°, b = —89.3°) data assuming the heavy disk. 

We use a power law IMF such as <f> oc M~( 1+x \ For the disk stars we use the Scalo's 
IMF (Scalo 1986 [37]), which is flat, i.e., x = -1 lor O.O8M < M < 1M @ and x = 2.3 
lor IMq < M < 5OAf . For spheroid stars we take Richer & Fahlman's value x = 3.5 lor 
M m i n < M < 1M & . For dark halo, x is a parameter and mass varies as M m i n < M < 1M & . 
We put the low-mass end ol IMF as M m i n = O.OIM©. This value is suggested Irom 
Sutherland (1994)'s estimate ol the mass ol MACHOs and a star formation theory ol 
Palla, Salpeter & Stahlar (1983) [32]. The observations ol MACHO (Alcock et al. 1995 
[5]) and EROS collaboration (Aubourg et al. 1995 [9]) suggest that the mass ol MACHOs 
is not so much less than 0.1M & . 

In Figure 8 we show N as a lunction ol f3 for given x for the direction ol Tyson(1988)'s 
SGP (/ = 312.5°, b = —89.3°) data. N is not sensitive to f3 but strongly depends on x. This 
is because N is determined mainly by the number ol hydrogen burning limit stars in the 
halo at the distance ol D = 3 ~ 5kpc. In this distance, the structure ol the halo changes 
little as f3 changes, as discussed in section 3. The number ol the shining low-mass stars 
is very sensitive to x. For other Tyson(1988)'s fields the figure does not changes so much. 
Since Tyson (1988) observed about 50 stars in a field, we conclude that IMF's power in 
the halo is steep at least as x ~ 5. 
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Recently, using HST, Bahcall et al. (1994) [12] observed a field of 4.4arcmin 2 wide, with 
I < 25.3 and V < 25.6 - 0.3(V - /). They found five stars with 2 < V — I (M < O.4M 
if brown dwarfs). Size of their field is about half of Tyson's one. Similar calculations show 
x > 6. 

5 Conclusion 

We investigated how Galactic mass models and the flattening of the halo affects r under 
the assumption that the halo consists of only MACHOs. We found that the constraints 
of the inner rotation curve almost uniquely determines the structure of the region of the 
halo (< 20kpc) where the microlensing events mainly occur. As a result r for LMC 
events varies at most a factor 2.5 from the standard spherical flat rotation curve model of 
r ~4xl0~ 7 . For / ~ 1, i.e., the halo consists of only MACHOs, the observation of MACHO 
collaboration 8tg 4 X 10~ 8 is consistent only if we take the most centrally concentrated model 
with f3 ~ 0.6, which means the rotation curve at r ~ 20kpc is declining. 

We also showed that /3-dependence of tsmc I tlmc is large for e > 0.3. This implies that 
tsmc I tlmc does not necessarily determine the flattening of the halo unless the rotation 
curve at r ~ 20kpc is determined definitely. 

For the bulge events, r cannot reach the observational value even in our models in 
which effects of the ellipticity and non flat rotation curve are included. This means that 
we must consider other effects such as a bar in the bulge. As Gould(1994a) [20] suggested, 
all of the matter inside the Sun may be in the thin disk and the LMC event may be by 
the stars in LMC. Anyway we need further observation and statistics (e.g. Gould (1994b) 
[21], Miyamoto & Yoshii (1994) [29], Gould (1994c) [22], Ansari et al. (1995) [6], Gould 
(1994d) [23]). 

For LMC events, we calculated the power x of IMF of MACHO consistent with Tyson's 
CCD survey as well as Bahcall et al. 's observation by HST. It is found that x is greater 
than 5. This suggests that the halo IMF is essentially a ^-function if the halo consists of 
only low mass stars. High value of x is required even if / ~ 0.2 unless the microlens events 
are occurred near stars in LMC. 

We thank Professor H. Sato for useful discussions. This Work was supported by Grant- 
in- Aid of Scientific Research of the Ministry of Education. 
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